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(57) Abstract: A volume rendering process is disclosed for improving the vi- 
sual quality of images produced by rendering and displaying volumetric data in 
voxel format for the display of three-dimensional data on a two-dimensional dis- 
play, with opacity (1402) and shading (1404) to control the realistic display of 
images rendered from the voxels. The process includes partitioning the plurality 
of voxels among a plurality of slices with each slice corresponding to a respective 
region of the volume. Each voxel includes an opacity value adjusted by applying 
an opacity curve to the value. The opacity value of each voxel in each cell in the 
volume is converted into a new visual opacity value that is used to calculate a new 
visual opacity gradient (1403) for only one voxel in the center of each cell. The 
visual opacity gradient is used to calculate the shading used to modify the color 
of the individual voxels based on the orientation of opacity isosurfaces passing 
through each voxel in the volume in order to create a high quality realistic image. 
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METHOD AND APPARATUS FOR VISUALIZATION OF 3D 
VOXEL DATA USING LIT OPACITY VOLUMES WITH SHADING 


10 


Field of the Invention 

This invention relates generally to the field of computer graphics. Particularly, . 
1 5 this invention relates to volume rendering. More particularly, this invention relates to the 
display of three-dimensional (3D) data on a two-dimensional (2D) display with shading 
and opacity to control the realistic display of volumetric data in voxel format. 


20 


Background of the Invention 


Volume rendering is an important area of computer graphics. It is employed in a 
wide variety of disciplines, including medicine, geology, biology and meteorology. 
Volume rendering allows a user to look inside an object and see features that were 
otherwise shielded by the rendering of the surface features of the object. One patent 
25 teaching volume rendering using voxels is U.S. Patent No. 6304,266, issued October 16, 
2001, entitled "Method and Apparatus For Volume Rendering" which is incorporated 
herein by reference. 

In Fig. 1 volumetric data is shown as consisting of a three-dimensional (3D) 
30 dataset of elements called "voxels" 102. Typically, the voxels 102 are uniformly 

distributed throughout a volume 104. Each voxel 102 has a position in the volume and 
has associated with it information such as color, illumination, opacity, velocity, 
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amplitude, etc. The information associated with each voxel 102 is produced by such 
disciplines as medicine (e.g., CAT scans), biology (confocal microscopy), and geoscience 
(seismic data). 

5 In Fig. 2 is shown how data values of voxels 1 02 are typically stored in a storage 

array 202. The position of a particular voxel in the volume is inherent in its location in 
the array. For example, array position 204 might be associated with a point 106 (Fig. 1) 
in the volume that is a specified distance from a specified comer of the volume. 
Typically, a single value is stored in the array 202 for each voxel 102, although it is also 
1 0 possible to store more than one value for each voxel 1 02, such as for color, illumination, 
opacity, velocity and amplitude. 

Figs. 3, 4, and 5 show 2D texture rendering subdividing volumetric data into 
slices. 2D texture rendering organizes the slices into three sets of slices 302, 402 and 502 
1 5 along, three different orthogonal axes. The voxels are partitioned among the sets of slices 
302, 402 and 502 and into cells containing multiple voxels in each slice. The partitioning 
is done based on the position of the voxels in array 202. In a similar fashion, 3D texture 
rendering typically slices the volume perpendicular to the viewing direction. 

20 The rendering is then accomplished on a slice-by-slice basis, moving from the 

rear-most slice 304, 404 and 504, respectively, to the front-most slice 306, 406 and 506 
respectively. The set of slices that is chosen and processed is the set whose axis makes 
the smallest angle to the viewing direction. While a new image is rendered it is blended 
with the previously drawn scene, creating perception of a 3D body. 2D texture rendering 

25 organizes the slices along one of the three volume dimensions, while 3D texture 

rendering slices the volume perpendicular to the viewing direction, which improves 
image quality, but requires interpolation between the volume data points. Such 
interpolation is usually performed by specialized graphics hardware. 
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Figs. 6 and 7 show how texture value, or "texel," is determined for each voxel in 
each slice (blocks 702 and 704). The texels are stored in a data buffer 602 (block 706). 
Typically, the texel value is an indication of the RGB colors (red, green & blue) to be 
displayed for a voxel as determined by one or more parameters dependent on the data 
5 value or values associated with the voxel and is found in a look-up table. For example, 
the texel data may include a value for each of the red, green, and blue (RGB) components 
associated with the voxel. When all of the voxels in the slice have been processed (block 
704), the contents of the data buffer are downloaded into a texture memory 604 (block 
708). 

10 

Fig. 8 shows a display device 802 upon which information downloaded with the 
texel data is displayed. Based on that information and the perspective requested by the 
user, the display device maps the texels onto pixels on a display screen 804 (block 710). 
As each slice is downloaded and rendered, the user sees the volume in the requested 
15 perspective. Each time the user changes the view, for example by using a software tool 
to rotate, translate or magnify the image of the volume, the process of downloading and 
rendering slices is repeated. 

In Fig. 9 is illustrated the display of a volume that shows the outside surfaces of 
20 the volume. The interior of the volume is not seen. 

In some applications, greater flexibility is achieved by using semi-transparent 
data. Semi-transparent data is created by adding an additional factor, alpha (a), to each 
voxel along with the RGB (red, green & blue) components described above. The value 
25 of alpha of a voxel determines the opacity of the voxel. Opacity is a measure of the 

amount a particular voxel on a slice will allow a voxel on a background slice that maps to 
the same pixel on a 2D display to show through. The opacity of a voxel controls how the 
image of the voxel is blended with the images of the voxels behind it in the view being 
displayed. An opacity value of 0.0 means a voxel is completely transparent and cannot 
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be seen so has no effect on the color of the displayed pixel on the 2D display since it is 
considered to be empty; and a value of 1 .0 means the voxel is completely opaque, may be 
considered solid and, if it has no other voxels mapped in front of it, its texel determines 
the color of the displayed pixel. Intermediate opacity values correspond to intermediate 
levels of opacity, and the texel defined colors of two voxels mapped to the same pixel are 
mixed in conventional ways to determine the color of the pixel that will be displayed. 

In Fig 10 is illustrated an opacity tool, such as the one included in the GEO VIZ 
product from Schlumberger-GeoQuest, the assignee of the present invention. The 
opacity tool is used to map volumetric data, such as geophysical seismic interpretation, 
magnetic imaging, and ultrasonography data, to see semi-transparent volumetric data 
therein. In those cases, the value of each voxel is not only mapped to a color defined by 
its texel but also with an opacity defined by alpha (a). In Fig. 10, the user has adjusted 
the opacity mapping, shown graphically by curve 1002, to make transparent all voxels 
(alpha = 0) except those having large positive or negative values. This has the effect of 
making most of the data transparent when displayed, as can be seen from the histogram 
1004 that reflects the distribution of the values of the voxels in the volumetric data 
displayed in Fig. 9. This is a technique called "standard opacity volume rendering" and 
allows a user to make voxels within a selected range of data values invisible, while 
leaving other voxels visible. 


In Fig. 1 1 is shown a display that results when the data displayed in Fig. 9 is 
processed using the opacity tool shown in Fig. 10. The surface of the volume no longer 
obscures structures inside the volume as is evident when comparing Figs. 9 and 1 1 . 

It is also apparent from the histogram 1004 and Fig. 1 1 that most of the opacity- 
adjusted voxels are transparent and have no effect on the display. 
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Summary of the Invention 

The invention is concerned with improving the visual quality of images produced 
by rendering volumetric data in voxel format for the display of three-dimensional (3D) 
5 data on a two-dimensional (2D) display with shading and opacity to control the realistic 
display of images rendered from the voxels. 

One of the main shortcomings of standard opacity volume rendering of the voxel 
format data volumes described in the Background of the Invention is that objects 
1 0 displayed using a voxel data volume appear flat, which inhibits depth perception and 
makes it hard, if not impossible, to determine the 3D shape, orientation, and relative 
positions of the objects. 

Figure 12 shows the display of an elliptical object rendered using standard opacity 
1 5 volume rendering of volumetric data consisting of voxels. Although it is clear that the 
cross-sectional shape is oval, there is no way to tell if the object is elliptical, ovoid, 
conical, etc. It is also impossible to tell if the object touches the flat seismic section 
behind it. The reason for these shortcomings is that traditional computer graphics 
surface lighting techniques do not work for voxels because the individual voxels have no 
20 proper surfaces to be lit. 

Figure 13 shows how an elliptical set of voxels will appear when displayed as a 
set of solid "bricks" lit with a light. The actual shape of the rendered object can now be 
seen. 


25 


In an attempt to alleviate the shortcomings with rendering and displaying voxels, 
the viewer makes elaborate adjustments to the volume color map or the opacity values, 
watches the rendering progressing from back to front, moves the viewpoint or simply 
employs the imagination to help recreate the objects in his or her mind. The methods are 
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very time-consuming, subjective and difficult for an inexperienced viewer to master, and 
they all yield unsatisfactory results. 


5 one or more light sources. When rendered and displayed the lighted volume offers the 
viewer ample visual infonnation to aid in the perception of depth, as well as the shapes, 
orientations, and positions of objects in the voxel volume. Lighting parameters are 
computed, and graphical attributes of individual voxels are adjusted based on the 
orientation of opacity isosurfaces passing through each voxel in the volume. 


To improve the display of objects the opacity of displayed voxels of a data 
volume may be adjusted depending upon the opacity of nearby voxels. Using an opacity 
tool each voxel in a volume is assigned a standard opacity value a using a 3-dimentional 
opacity function a(ij,k) where alpha represents opacity and the letters ij,k represent 
15 orthogonal directions. 

Isosurfaces connect voxels with equal opacity values. For each voxel, the 
algorithm estimates a normal to the isosmface passing though the voxel, which is equal to 
the negative gradient of the opacity function at the voxel, and uses it to shade the voxel as 
20 if it were a point on the isosurface. Though results are typically shown for a single uni- 
directional white light source, the results can be easily extended to bi-directional light 
sources or colored lights. 

In Fig. 14 are shown the basic steps implemented in the rendering process that 
25 includes the teaching of the present invention, and these steps are generally described 

immediately below. The invention is implemented only in the second and third steps, in 
blocks 1402, 1403 and 1404. 
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At block 1401, the initial opacity value of each voxel in the input voxel data 
volume is processed with an opacity tool, such shown in Fig. 10, to make most of the data 
transparent when displayed, as previously described, and the structures inside the volume 
are no longer obscured as can be seen in Fig. 1 1. As part of this process the initial 
5 opacity value for each visible voxel in the data volume is converted to a standard opacity 
(a) value. 

At block 1402, the standard opacity value a for each visible voxel is converted to 
a new "visible opacity" p value in accordance with the teaching of the present invention. 

10 

At block 1403, a new visible gradient value is calculated for every visible voxel in 
a volume using the new visible opacity p values in accordance with the teaching of the - 
present invention. The new visible gradient value calculated for each voxel accounts for 
degenerate cases in numeric gradient computation and are used in all further rendering 
1 5 computations. Only visible cells, selected using the opacity tool and having a new visible 
opacity p calculated, have their opacity p value used to compute a new visible gradient, 
G , for each voxel. The new visible gradient, G , for each voxel are used to shade and 
render the displayed voxels. 

20 At block 1404, ambient and diffuse shading is computed that would be applied to 

an opacity isosurface passing through the voxels in an ambient and diffuse illumination 
model wherein the voxel volume is illuminated with one or more light sources (typically 
directional and bi-directional). The direction of the negative visible gradient vector 
G serves as a normal to the isosurface. Special treatment is added for voxels inside 

25 opaque areas of the volume based on specifics of volumetric geoscience data, which 

improves a user's perception of the rendered image. A lighted and shaded volume offers 
the viewer ample visual information to aid in the perception of depth, as well as the 
shapes, orientations, and positions of objects in the volume. 
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At block 1405, for both 2D and 3D textures, the colors of all visible voxels in the 
volume are modified by applying the shading computed in block 1404. 

At block 1406 the processed, visible, voxel data volume, with shading, is 
5 displayed. 

In one aspect of the invention, a computer is used to process voxel volume data, 
and the invention comprises an article of manufacture comprising a medium that is 
readable by the computer and the medium carries instructions for the computer to 
1 0 perform the novel process of calculating the new visible opacity and the visible gradient 
for each voxel as described in the previous paragraphs. 

Brief Description of the Drawings 

1 5 The invention will be better understood upon reading the following Description of 

the Preferred Embodiment in conjunction with the drawings in which: 
Fig. 1 is a representation of voxels within a cubic volume; 
Fig. 2 is a representation of the storage of voxel values in a memory; 
Figs. 3-5 illustrate the partitioning of a volume into slices in orthogonal planes for 
20 2D texture rendering; 

Fig. 6 is a representation of a volume rendered image without shading; 

Fig. 7 is a block diagram of the steps in a volume rendering process; 

Fig. 8 is a block diagram of volume rendering equipment; 

Fig. 9 is a rendering of an opaque volume; 
25 Fig. 1 0 is a representation of an opacity tool; 

Fig. 1 1 is a rendering of a semi-transparent volume after an opaque volume is 
processed using an opacity tool; 

Fig. 12 shows the display of an elliptical object rendered using standard opacity 
volume rendering of volumetric data consisting of voxels without shading; 
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Fig 13 shows how an elliptical set of voxels would appear when displayed as a set 
of solid "bricks" lit with a light; 

Fig. 14 shows in block diagram form the steps of processing and displaying 
volumetric geoscience data consisting of a volume of voxel data; 

Figs. 15 - 17 are exploded views of volumes each showing a cell made up of 
voxels that is used to calculate new visual opacity and visual gradient values for voxels; 

Figs. 1 8- 20 are tables showing how the new visual opacity values for voxels are 
calculated; and 

Figs. 21 - 23 are tables showing how the new visual opacity values are used to 
calculate new visual gradient values for voxels. 

Description of the Preferred Embodiment 

In the following description all reference numbers indicate in which Figure they 
15 are located. Thus, for example, reference number 1004 is found in Fig. 10, and reference 
number 1405 is found in Fig. 14. 

Volumetric data volume is first partitioned into slices as illustrated in Figs. 3, 4 
and 5, and the slices each contain rows and columns of voxels, not seen in these Figures 
20 but seen in Figures 15-17. Volume rendering of the data may be speeded by reducing 
the number of voxels that are processed, downloaded and rendered by eliminating 
transparent voxels. Such volume rendering is taught in the previously mentioned U.S. 
Patent No. 6,304,266 that is incorporated herein by reference. 

25 In Figs. 15 - 17 are shown a representative volume of voxels formed into a 

plurality of slices A - E and cells, as particularly applied to geoscience data, voxels are 
grouped into cells containing no fewer than twenty-seven alpha-numeric designated 
voxels (eg. Al 1, C21) and forming a cube. Other undesignated voxels shown in these 
Figures are other voxels of the data volume in which the cells being described are 
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located. Only one cell is shown in each of these Figures for simplicity, but in reality 
there would be many slices and many cells in each slice of a volume. These cells are 
described in greater detail hereinafter. The slices are shown separated from each other 
only for ease in seeing the twenty-seven voxels Ai i-A33, Bl 1-B33 and CI 1-C33 that 
5 make up cells one and three in Figs. 16 and 17 respectively, and the twenty-seven voxels 
Bl 1-B33, CI 1-C33 and Dl 1-D33 that make up cell two in Fig. 16. 

For the purposes of this example, only horizontal slices are specifically shown in 
Figs. 15, 16 and 17 but there are also vertical slices, and the intersecting orthogonal slices 

10 through the data volume are used to create the three dimensional rows and columns of 
voxels, as shown, so that, in this particular case, the vector / arrow direction of the light 
source is parallel to one of the three orthogonal axes («,v,m>) of the sliced data volume 
while the other two axes are orthogonal to the direction of the light source vector. In 
Figs. 15 and 16 the light vector is parallel to axes W, and in Fig. 17 the light vector is 

15 parallel to axes V. Note that while the illustrations show the light direction parallel to 
one of the axes, the .method employed supports arbitrary directions for the light source. 
The value of initial opacity, a, for the voxels in each cell are used to calculate the values 
of visible opacity p only for the voxels in each cell as described hereinafter with 
reference to Figs. 1 8 - 20. 

20 

While the preferred cell size, as shown in Figs. 15 - 17, is twenty-seven voxels 
forming a cube of 3x3x3 voxels, a different cell partitioning size might be appropriate for 
other data, such as medical or meteorological data. An example of another cell size may 
be one-hundred, twenty-five voxels forming a cube of 5x5x5 voxels. The number of 
25 voxels in a cell is based on the distribution of data in the volume. 

In Fig. 1 5 cell one the particular voxel of interest is the center voxel B22. In Fig. 
15 light is incident on the top of the data volume and cell one and is represented by the 
light source vector / arrow. In Fig. 16 cell two the particular voxel of interest is the center 
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voxel C22. In Fig. 16 light is incident on the top of the data volume and cell two and is 
represented by the light source vector / arrow. In Fig. 17 cell three is made up of the 
same voxels as cell one, but light is incident on the right side of the data volume and cell 
three as represented by the light source vector / arrow. 

5 

As briefly described above, the initial opacity value for each voxel in a data 
volume is first mapped with an opacity tool, such as the one illustrated in Fig. 10 and that 
is included in the GEO VIZ product of Schlumberger-GeoQuest, to obtain a "standard 
opacity" value a for each voxel in a data volume. The opacity tool adjusts the opacity 

10 mapping of voxels in the data volume to make transparent chosen voxels (a = 0) except, 
for example, those having large positive or negative values. This has the effect of 
making most of the data transparent when rendered, and structures inside the data volum^ , ; 
are no longer obscured as can be seen in Fig. 1 1 . This is a technique called "standard 
opacity volume rendering" and allows a user to make voxels within a selected range of 

15 data values invisible, while leaving others visible. The operator does this by changing the 
setting of the opacity tool. This technique is described in greater detail in U.S. Patent No. 
6,304,266 cited above. 

In a simplified treatment of visibility we, first, find the largest component of the 
20 light incident on a data volume and mark it as the "visibility direction" (e.g., «, assuming 
that L u = max(z, w ,Z v ,JL lv )). This would be the vector / arrow in Figs. 15-17. Thus, for 
a light coming from the u direction (not represented in Figs. 15-17), the new "visible 
opacity" p value is computed using the following equations: 
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P{i-lJ,k) = a{i-lJ,k) 
P(i, j , /r) = max(/?(i - i, y , A), y, A)) 
+ i,y,/c) = max(/?(/,y\ *),a(i + l,y\*)) 

+ 7,y,/c ) = a(i + i,y\/r) 
j3(i,j,k) = m*x(j3(i + lJ 9 kMiJ 9 kj) 
p{i - i, y , /c ) = max(^(i, y , A ), a(i-lj, k)) 

The above mathematical equations are used to calculate the visible opacity p for 
each voxel in each of the nine columns of three voxels in a cell, as viewed from the 
5 direction of the light, such as voxel column All, Bl 1, CI 1 and voxel column A23, B23, 
C23 in cell one. When a voxel column is viewed from the direction of the light source, if 
the first voxel is opaque all other voxels behind it cannot be seen. For example, in cell 
one in Fig. 15, if voxel Al 1 is opaque, voxels Bl 1 and CI 1 behind it cannot be seen. 
Similarly, if voxel Al 1 is semi-transparent (initial opacity value between 0.0 and 1.0) the 
10 opacity value of voxels Bl 1 and CI 1 behind voxel Al 1 cannot be any less. For example, 
if voxel Al 1 is semi-transparent with an opacity value of 0.7, voxels Bl 1 and CI 1 cannot 
have a visible opacity value any lower than 0.7 and are changed accordingly in 
calculating visible opacity values p for voxels Bl 1 and CI 1. 

15 Figs. 18-20 are tables used to simplify and more clearly show the mathematical 

calculations performed by the above equations to derive the new visible opacity values P, 
respectively, for all twenty-seven voxels in each of cells one, two and three using the 
"standard opacity" value a of the voxels in accordance with the teaching of the present 
invention. In the following description, "a" indicates standard opacity and when used as 

20 ctAl 1 indicates the standard opacity of voxel A 1 1, and "p" indicates the new visible 
opacity and when used as pAl 1 indicates the visible opacity of voxel Al 1. 
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Returning to Fig. 1 5 to describe calculating visible opacity p values for all voxels 
in cell one as shown in the table in Fig. 1 8. With light being incident on the top of the 
volume in Fig. 1 5, and the Al 1-A33 voxels being on the top of the cell, the standard 
opacity value of each of the voxels Al 1-A33 (aAl 1 - aA33) is equal to the visible 
5 opacity (3 value of each of these voxels. There is nothing in front of these voxels to block 
or alter the intensity of the light shining upon them, so their standard opacity will equal 
their visible opacity. This is represented in the "A slice" column in Fig. 1 8 as 


standard opacity of voxel Al 1 (aAl 1) equal to the visible opacity pAl 1 of voxel All, 
when p A 1 1 is greater than the standard opacity of voxel B 1 1 (aB 11) behind it, then the 
value of visible opacity of voxel Bl 1 (PB1 1) cannot be any smaller and is changed to 
equal the visible opacity of voxel Al 1. That is, the visible opacity of voxel Bl 1, pBl 1, is 

15 set equal to the visible opacity value of voxel Al 1, PAl 1 (pBl 1 = pAl 1). Conversely, if 
the visible opacity pAl 1 is less than or equal to the standard opacity aBl 1, then the 
visible opacity of voxel Bl 1, pBl 1, is set equal to its standard opacity value (pBl 1 = 
aBl 1). Continuing down the same column of voxels Al 1, Bl 1 and CI 1 to the CI 1 
voxel, and again using the same rationale, if the visible opacity of voxel Bl 1, pBl 1, is 

20 greater than the standard opacity of voxel CI 1 aCl 1, then the visible opacity of voxel 
Cll, PC11, is set equal to the visible opacity of voxelBll, pBll (pCll = PB11). 
Conversely, if the visible opacity of voxel Bl 1 5 pBl 1, is less than or equal to the standard 
opacity of voxel CI 1, aCl 1, then the visible opacity pCl 1 remains equal to its standard 
opacity (PCll =aCll). 


This same determination of visible opacity for all voxels in cell one is repeated for 
each of the other of the nine columns of three voxels in Fig. 15 (e.g. voxel columns A32, 
B32 & C32; A22, B22 & C22; etc.). The visible opacity p values calculated in this 


aAl l=pAl 1, aA12=pA12 etc. through aA33=pA33. 


10 


To calculate the visible opacity p of voxel Bl 1 behind voxel All, and with 


25 
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manner for all voxels in cell one are used for the sole purpose of calculating the visual 
opacity gradient for only voxel B22 in the center of cell one as is described hereinafter 
with reference to Fig. 21. Although other cells, such as cell two, include many of the 

Ramp w/*w<*lc J*-»^1n/-I^/-J '* ~ 1 „ . jl 1 . r~ , .. . 

— „~ , **,w*„«w« ui wiic, uic vaiucs ui visual opacity will be recalculated for all 
5 voxels in each cell, and the visual- opacity values for voxels in cell two most likely may 
be different than those values calculated when the voxels are in cell one. 

Although only one cell is shown in Fig. 15 there are in reality many cells in a data 
volume and the visible opacity p is determined for all voxels in the data volume. This is 
1 0 done by having each voxel at the center of a cell and performing the calculations 
described above. The voxels on the outer surface of the volume can be disregarded. 

In Fig. 19 is a table showing how standard opacity a values of each of the voxels 
in cell two (Fig. 16) are used to obtain the new visible opacity p values for each of the 
1 5 voxels in cell two. With light being incident on the top of the data volume in Fig. 1 6, 
and the Bl 1-B33 voxels of cell two being on the side of the cell from which the light is 
coming, the standard opacity value of each of the voxels Bl 1-B33 is equal to the visible 
opacity value of each of these same voxels. This is represented in the "B slice" column 
in Fig. 19 as aBl l=pBl 1, ctB12=pB12 etc. through aB33=pB33. 

20 

When calculating the visible opacity of voxel C 1 1 behind voxel B 1 1 , and using 
the rationale described in previous paragraphs, if the visible opacity pBl 1 is greater than 
the standard opacity a of voxel CI 1, then the visible opacity of voxel CI 1 is changed to 
equal the visible opacity of voxel B 1 1 . That is pCl 1 = pB 1 1 . Conversely, if the visible 
25 opacity of voxel B 1 1 , PB 1 1 , is less than or equal to the standard opacity of voxel C 1 1 , 
aCl 1, then the visible opacity pCl 1 of voxel CI 1 is set equal to its standard opacity. 
That is pCl 1 = aCl 1 . Continuing down the column of voxels Bl 1, CI 1 and Dl 1, if the 
visible opacity of voxel CI 1, pCl 1, is greater than the standard opacity of voxel Dl 1 
oDl 1, then the visible opacity of voxel Dl 1 , 0D1 1, is set equal to the visible opacity of 
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voxel Cll, pCll. ThatispDll = PC11. Conversely, if the visible opacity of voxel 
Cl l,pCl 1, is less than or equal to the standard opacity of voxel Dl 1, aDl 1, then the 
visible opacity PD1 1 remains equal to its standard opacity aDl 1. That is PD1 1 = aDl 1. 

5 This same calculation of visible opacity p for all voxels in cell two is repeated for 

each of the other nine columns of three voxels in cell two in Fig. 19 (e.g. voxel columns 
B32, C32 & D32; B22, C22 & D22; etc.). In this manner the visible opacity p of all 
voxels in cell two is calculated. The calculated value of visible opacity p for all voxels in 
cell two are only used in the equations in Fig. 22 to calculate the visible opacity gradient 
10 p of voxel C22 in the center of cell two. 

Although only one cell is shown in Fig. 16 there are in reality many cells in a data 
volume and the visible opacity p. is determined for all voxels in the data volume. This is 
done by having each voxel at the center of a cell and performing the calculations 
15 described above. The voxels on the outer surface of the volume can be disregarded. 

In Fig. 17 is cell three that has the same voxels as cell one shown as in Fig. 15 but 
the direction of light on the data volume and cell three is from the right side, rather than 
from the top. Accordingly, calculations for visible opacity p are identical in manner but 

20 are different. The columns of three voxels are on their side through cell three. For 

example, voxels B23, B22, B21 and voxels C13, C12, Cl 1. The tables for calculating 
visible opacity P for all voxels in cell two are shown in Fig. 20. In view of the previous 
description of how these calculations are made with reference to the tables in Figs. 18 and 
19, the description is not repeated here for the sake of brevity. The visible opacity p of 

25 all the voxels in cell three are used to calculate the visible opacity gradient G only for 
voxel B22 in the center of cell three as described hereinafter with reference to Fig. 23. 
Again, there are many cells in the data volume in Fig. 1 7 and visible opacity is 
determined for all voxels. 
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Figs. 21 - 23 are tables used to describe the mathematical calculations performed 
to derive the three gradient components G„ , G v , and G w that define a new visible 
opacity gradient G for only the voxel in the center of each of representative cells one, 
two and three in accordance with the teaching of the present invention. The gradients 
must be derived for all voxels in a data volume so similar tables are derived for the other 
voxels in a data volume, not just cells one, two and three. The three gradient components 
are calculated for every voxel using the newly calculated value of visible opacity p for all 
voxels in each cell. For cell one in Fig. 15 the center voxel is B22; for cell two in Fig. 16 
the center voxel is C22; and for cell three in Fig. 17 the center voxel is B22. The new 
visible opacity gradient G for all voxels are then used to render the voxel data volume in 
a manner well known in the art. These calculations are repeated for each voxel in a - 
volume and are described in greater detail hereinafter. 

Returning to cell one to describe how the visible opacity p values for each of the 
1 5 twenty-seven voxels in cell one are used to calculate the new visible opacity gradient G 
for center voxel B22. Using the visible opacity p values calculated for each of the 
twenty-seven voxels Al 1 through C33 in cell one, these values are used in the equations 
shown in Fig. 1 8 to calculate vector components G u , G v and G w of vector G only for 
voxel B22 in the center of cell one. The vector components are then combined to get 
20 vector G , the negative visible opacity gradient for voxel B22. 

In the same manner, the new visible opacity P values calculated for each of the 
twenty-seven voxels Bl 1 through D33 in cell two are used in the equations in Fig. 19 to 
calculate gradient components G u , G v and G w of the visible opacity gradient vector 
25 G only for voxel B22 in the middle of cell two. The gradient components are then 
combined to get the negative visible opacity gradient for voxel B22. 


5 
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Also, the visible opacity p values for the twenty-seven voxels Al 1 through C33 in 
cell three are used in the equations in Fig. 20 to calculate vector components 
G u , G v and G w of the visible opacity gradient vector G only for voxel C22. The gradient 
. components are then. combined to get the negative visible opacity gradient for voxel C22. 
5 The basic gradient equations are known in the prior pail and are described in more detail 
in the following paragraphs, but values of visible opacity p are used in the equations 
rather than values of standard opacity as previously used to derive the simplified 
calculations shown in Figs. 21-23. 

10 A negative opacity gradient, G , at a particular voxel is determined by three 

partial derivatives along the three major axes, G u , G v , and G w as: 

G =(- G u ~ G v i-G w ) 
where the standard gradient equations used are: 

G v = y Za(i,j + l,k)- y £ i a(iJ-l,k) 

i,k i,f( 

G w ^a(i,j,k + l)-Y J oc(i,j,k-l) 

iJ iJ 

1 5 However, in accordance with the teaching of the invention, these three standard 

gradient equations for calculating vector components G ti , G v andG^ are modified to 

calculate a new negative "visible opacity gradient" vector G , shown in the following the 
equations, by using the new visible opacity p values, rather than standard opacity values 
a shown in the equations immediately above. Substituting a value of visible opacity p in 
20 the standard gradient equations, the gradient equations then become: 

G v =YP( i >J + 1 > k )-llP( i >J- 1 > k ) 
'J i,j 
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where p(ij,k) are the visual opacity of a single voxel along each its three major axes ij,k. 

Only visible voxels, selected using the opacity tool, and then further processed to 
derive visible opacity p for each visible voxel, are used lo compute vector components 
G u , G v andG w of negative visible opacity gradient G , preferably using the 26- 
neighborhood central difference method described above, at each visible voxel using the 
modified equations immediately above. The negative visible opacity gradient G is 
calculated using the equation: 

where the vector components G u , G v andG w are calculated using the tables in Figs. 21 - 
23. These tables reflect the calculations in the above gradient equations. 

Gradient G computed using the twenty-six neighborhood difference method have 
significantly more distinct values and result in smoother images than a six-neighborhood 
difference, although the latter is faster to compute but has only twenty-seven distinct 
values that result in abrupt changes of shading. For the present description we use the 
twenty-six neighborhood difference gradient of the opacity. These equations work 
regardless of the number of lights illuminating an object or the color of the light. 

The new visible opacity gradient, G , for each of voxels B22 (cell one), C22 (cell 
two) and B22 (cell three) are then used to calculate ambient and diffuse shading intensity 
for those voxels in a manner previously known, with the addition of special treatments . 
provided by this invention, and generally described with reference to block 1404 in Fig. 
14. Shading is computed as if it were applied to an opacity isosurface passing through all 
25 voxels in a volume in an ambient and diffuse illumination model wherein the voxel 
volume is illuminated with one or more light sources (typically directional and bi- 
directional). The direction of the negative visible opacity gradient vector serves in this 
case as a normal to the isosurface. Special treatment is added for the voxels inside 
opaque areas of the volume based on specifics of volumetric geoscience data, which 
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improves a user's perception of the rendered image. A lighted and shaded volume offers 
the viewer ample visual information to aid in the perception of depth, as well as the 
shapes, orientations, and positions of objects in the volume. 

5 As part of computing shading a decision is made if the volume being rendered is 

to be shaded as if lit by a uni-directional light source pointing in one direction or by bi- 
directional light consisting of two identical directional lights pointing in opposite 
directions on a volume. When the volume being rendered is lit by a uni-directional 
source the ambient and diffuse shading intensity is calculated using opacity 
10 (7^=/(^)) as follows: 

ifGL>0 

1 shotting = 1 ambient + 1 diffuse i G ' L )> Norm(G ) 

else if G L<0 

^ shading ~~ ^ ambient 

where G is the negative visible opacity gradient, calculated as previously described, and 
L is the unit light vector. 

For bi-directional light consisting of two identical directional lights pointing in 
opposite directions on the volume, the ambient and diffuse shading intensity is calculated 
as follows: 

ifGL>0 

I shading = 1 ambient + I diffuse ( G ' L ) l Norm(G ) 

else if G L <0 

1 shading = 1 ambient ~ I diffuse (p ' I*) 1 Nomt{G ) 

where G is the vector of the negative visible opacity gradient and L is the unit light 
vector. 
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For the case of G = 0 (when and only when Norm(G) = 0 ) deserves special 
treatment. There exist four possible scenarios: 

a) G = 0 and the opacity value in the voxel is 0.0 ( J3{i, j, /c ) = 0 ). We treat this 
5 case as an empty space and set I shadiItg = 0 . 

b) G - 0 while the data value in the voxel is not 0.0 (J3(i,j 9 k) data values 
of the surrounding voxels are not all 0.0, and 


10 


jjk l,k 


In this case the voxel is treated as positioned within a homogeneous voxel body. 
Theoretically, in this case an isosurface that would pass through the voxel is not defined 
and, thus, the voxel should not be shaded ( I s/ia(Jing = 1 ). This produces visual artifacts 
when somebody observes the rendered image. In order to eliminate it, we choose to 
1 5 assign such voxels an artificial gradient G prefered and set the shading intensity in such a 
voxel to be: 

I shading ~ ^ ambient ^ ? diffuse (gprefered * ^ ) 

For most of 3D geoscience volumes there are horizontal layers in the data so the 

is: 

G prefered = {o,o a) 


20 most reasonable choice for G prefered is: 


c) G =0 while the data value in the voxel is not 0.0 (J3(ij,k) * 0 ), but data 
values of all the surrounding voxels are 0.0. In this case of a single voxel surrounded by 
25 empty space the direction of the visual opacity gradient in not defined. Thus, we 
arbitrary select it to satisfy the equations from scenario (b) immediately above. 
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d) The rest of the scenarios might require re-computing G by using right 
differences or some other known method. The required computations are time 
consuming. Hie shading computed using the formula in scenario (b) above for most 
5 voxels approximates the precisely computed ones reasonably well. 

After the shading intensity has been computed, including the special treatment in 
scenarios (a) - (d), it is applied to a color-coded opacity volume. This was referred to 
with reference to block 1405 in Fig. 14. There are multiple ways known to implement 
10 this. Often, the color of a voxel is derived from the associated data value, // , using a 
color palette, also called a look-up table or LUT, as: 

(r,g,b 9 a)=LUT(ju) 

where LUT is essentially a one-dimensional array of r,g,b,a quadruples indexed by the 
data value, fi . Thus, whenever a color palette is used, the color of a voxel is a function 
15 of only the data value. 

When we apply shading to a voxel, the final color of the voxel depends not only 
on the data value, but also on the shading intensity associated with the voxel. The initial 
formula calls for multiplication of each of the color components of each voxel by the 

20 shading intensity for the voxel. This requires many additional computations, slowing the 
process. One could alternatively use an extended color palette, SHADINGJLUT, which 
is a two dimensional matrix composed of r,g,b,a columns computed for different values 
of shading intensity, I shading . Once pre-computed, such an extended color palette could 
be used to look-up a color of a shaded voxel using the data value, // , and shading 

25 intensity, I shading , as indexes: 

(r,*Aa)= SHADING _LUT(juJ sha(ting ) 
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The initial coloring (r,g,b) of a 3D seismic opacity volume does not depend on 
the volume opacity and can be specified based on any data values associated with a 
volume (e.g., amplitude, instantaneous frequency, etc.). In order to simulate illumination 
of a voxel by a white light, the color components of voxel material are modified by 
5 multiplying them with the shading intensity: 

(''shaded » S shaded > b 'shaded >°0 = C shading ■' »^ shading §■> 1 shading A«) 

The results obtained from this calculation are the rendered volumetric data that is used to 
display the 3D seismic information volume on a 2D display device. 

1 0 While what has been described hereinabove is the preferred embodiment of the 

invention it will be appreciated by those skilled in the art that numerous changes may be 
made without departing from the spirit and scope of the invention. 


*f* H* H* 


What is claimed is: 
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CLAIMS 

1 1. A method for rendering a volume of voxel data with shading and opacity, wherein 

2 each voxel comprises a value representative of a parameter at a location within the 

3 volume and each voxel has an initial value of opacity, the method comprising the steps 

4 of: 

5 calculating a revised value of opacity for each of the voxels in the volume, the 

6 revised value of opacity of a voxel being dependent upon its initial value of opacity and 

7 the revised value of opacity of voxels proximate thereto; and 

8 calculating an opacity gradient for each of the voxels in the volume using the 

9 calculated revised values of opacity for all voxels. 


1 2. The volume rendering method of claim 1 wherein the step of calculating a revised 

2 value of opacity further comprises the step of: 

3 creating a cell of voxels surrounding each voxel in the volume, wherein all voxels 

4 in each cell are arranged into groups so that each voxel in each group of voxels within 

5 each cell are positioned one behind the other in a line parallel to the primary direction of 

6 light. 

1 3. The volume rendering method of claim 2 wherein the step of calculating a revised 

2 value of opacity further comprises the steps of: 

3 setting the revised value of opacity of the voxel closest to the source of light in 

4 each group of voxels in each cell equal to its initial value of opacity; and 

5 setting the revised value of opacity of all other voxels in each group of voxels in 

6 each cell equal to the revised value of opacity of an adjacent voxel in the same group of 

7 voxels that is closer to the source of light if the revised value of opacity of the closer 

8 voxel is equal to or greater than the initial value of opacity of the adjacent other voxel, 

9 and setting the revised value of opacity of the adjacent other voxel equal to its initial 
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10 value of opacity if the revised value of opacity of the closer voxel is less than the initial 

1 1 value of opacity of the adjacent other voxel. 

A • • iiiv vuiumv/ lwiiuvmig iii^inwi WX viaiiii VV 1J.OJL tnw V/A v^aiouKamig, till 

2 opacity gradient for each of the voxels in the volixme using the calculated revised values 

3 of opacity for all voxels further comprises the steps of: 

4 combining the revised values of opacity for the voxels in each cell to derive three 

5 orthogonal opacity gradient components for the voxel in the center of each cell; and 

6 combining the three orthogonal opacity gradient components for the voxel in the 

7 center of each cell to derive an opacity gradient that is normal to an isosurface passing 

8 through the voxel in the center of each cell. 

1 5. The volume rendering method of claim 1 wherein the step of calculating an 

2 opacity gradient for each of the voxels in the volume using the calculated revised values 

3 of opacity for all voxels further comprises the steps of: 

4 combining the revised values of opacity for the voxels in each cell to derive three 

5 orthogonal opacity gradient components for the voxel in the center of each cell; and 

6 combining the three orthogonal opacity gradient components for the voxel in the 

7 center of each cell to derive an opacity gradient that is normal to an isosurface passing 

8 through the voxel in the center of each cell. 

1 6. The volume rendering method of claim 4 further comprising the steps of: 

2 calculating shading for the volume using the opacity gradient; and 

3 displaying the rendered volume on the display device. 

1 7. The volume rendering method of claim 1 further comprising the steps of: 

2 calculating shading for the volume using the opacity gradient; and 

3 displaying the rendered volume on the display device. 
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1 8. A method for rendering a volume of voxel data with shading and opacity 

2 dependent upon the direction of a source of light specified to be illuminating the volume 

3 for rendering, wherein each voxel comprises a value representative of a parameter at a 

4 location within the volume and each voxel has an initial value of opacity, the method 

5 comprising the steps of: 

6 (1) calculating a revised value of opacity for a first voxel in the volume by creating a 

7 cell of voxels surrounding the first voxel, wherein all voxels in the cell are arranged into 

8 groups so that each voxel in each group of voxels within the cell are positioned one 

9 behind the other in a line parallel to the primary direction of light, the revised opacity 

1 0 calculation step further comprising the steps of: 

11 (a) setting the revised value of opacity of the voxel closest to the source of 

12 light in each group of voxels in the cell equal to its initial value of 

1 3 opacity; and 

14 (b) setting the revised value of opacity of other voxels in the groups of voxels 

1 5 in the cell equal to the revised value of opacity of an adjacent voxel in the 

16 same group of voxels that is closer to the source of light if the revised 

1 7 value of opacity of the closer voxel is equal to or greater than the initial 

1 8 value of opacity of the other voxel, and setting the revised value of 

1 9 opacity of the other voxel equal to the initial value of opacity of the other 

20 voxel if the revised value of opacity of the closer voxel is less than the 

21 initial value of opacity of the other voxel; 

22 (2) repeating the revised opacity calculations step for all other voxels in the volume 

23 other than the first voxel; and 

24 (3) calculating a gradient of opacity for each of the voxels in the volume using the 

25 calculated revised values of opacity for all voxels. 

1 9. The volume rendering method of claim 8 wherein the step of calculating an 

2 opacity gradient for each of the voxels in the volume using the calculated revised values 

3 of opacity for all voxels further comprises the steps of: 
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combining the revised values of opacity for the voxels in each cell to derive three 
orthogonal opacity gradient components for the voxel in the center of the last mentioned 
cell; and 

P.omVnrnncr thp> tVnv**» nrtVir\rrn»-ial rvr\o/-»i+ir m*r\Ai r^^-t- ~ ~ — -C — i i__ _ 
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center of each cell to derive an opacity gradient that is normal to an isosurface passing 
through the voxel in the center of the last mentioned cell. 


The volume rendering method of claim 8 further comprising the steps of: 
calculating shading for the volume using the opacity gradient; and 
displaying the rendered volume on the display device. 


1 1 . The volume rendering method of claim 9 further comprising the steps of: 
calculating shading for the volume using the opacity gradient; and 
displaying the rendered volume on the display device. 

12. A method for rendering a volume of voxel data with shading and opacity 
dependent upon the direction of a source of light specified to be illuminating the volume 
for rendering, wherein each voxel comprises a value representative of a parameter at a 
location within the volume and each voxel has an initial value of opacity, the .method 
comprising the steps of: 

(1) calculating a revised value of opacity for a first voxel in the volume by 
creating a cell of voxels surrounding the first voxel, wherein all voxels in the cell are 
arranged into , groups so that each voxel in each group of voxels within the cell are 
positioned one behind the other in a line parallel to the primary direction of light, the 
revised opacity calculation step further comprising the steps of: 

(a) setting the revised value of opacity of the voxel closest to the source 
of light in each group of voxels in the cell equal to its initial value of opacity; 

(b) setting the revised value of opacity of other voxels in the groups of 
voxels in the cell equal to the revised value of opacity of an adjacent voxel in the 


26 


WO 02/50779 


PCT/US01/48439 


15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 

1 

2 
3 

1 

2 
3 
4 
5 
6 
7 


same group of voxels that is closer to the source of light if the revised value of 
opacity of the closer voxel is equal to or greater than the initial value of opacity of 
the other voxel, and setting the revised value of opacity of the other voxel equal to 
the initial value of opacity of the other voxel if the revised value of opacity of the 
closer voxel is less than the initial value of opacity of the other voxel; 

(2) repeating the revised opacity calculations step for all other voxels in the 
volume other than the first voxel; 

(3) calculating three orthogonal opacity gradient components for the first voxel, 
one of the orthogonal opacity gradient components being parallel to the primary direction 
of light, the opacity gradient calculation step further comprising the steps of: 

(a) combining the revised values of opacity for the voxels in the cell to 
derive the three orthogonal opacity gradient components for the first voxel; and 

(b) combining the three orthogonal opacity gradient components for the 
first voxel to derive an opacity gradient that is normal to an isosurface passing 
through the first voxel; 

(c) repeating the opacity gradient calculation step for all other voxels in 
the volume other than the first voxel. 

13. The volume rendering method of claim 12 further comprising the steps of: 
calculating shading for the volume using the opacity gradient; and 
displaying the rendered volume on the display device. 

14 Apparatus for rendering a volume of voxel data with shading and opacity, 
wherein each voxel comprises a value representative of a parameter at a location within 
the volume and each voxel has an initial value of opacity, the rendering apparatus 
comprising: 

means for calculating a revised value of opacity for each of the voxels in the 
volume, the revised value of opacity of a voxel being dependent upon its initial value of 
opacity and the revised value of opacity of voxels proximate thereto; and 
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8 means for calculating an opacity gradient for each of the voxels in the volume 

9 using the calculated revised values of opacity for all voxels. 

1 15. The volume rendering apparatus of claim 14 wherein the means for calculating a 

2 revised value of opacity further comprises: 

3 means for creating a cell of voxels surrounding each voxel in the volume, wherein 

4 all voxels in each cell are arranged into groups so that each voxel in each group of voxels 

5 within each cell are positioned one behind the other in a line parallel to the primary 

6 direction of light. 

1 16. The volume rendering apparatus of claim 1 5 wherein the means for calculating a 

2 revised value of opacity further comprises: 

3 means for setting the revised value of opacity of the voxel closest to the source of 

4 light in each group of voxels in each cell equal to its initial value of opacity; and 

5 means for setting the revised value of opacity of all other voxels in each group of 

6 voxels in each cell equal to the revised value of opacity of an adjacent voxel in the same 

7 group of voxels that is closer to the source of light if the revised value of opacity of the 

8 closer voxel is equal to or greater than the initial value of opacity of the adjacent other 

9 voxel, and setting the revised value of opacity of the adjacent other voxel equal to its 

10 initial value of opacity if the revised value of opacity of the closer voxel is less than the 

1 1 initial value of opacity of the adjacent other voxel. 

1 1 7. The volume rendering apparatus of claim 1 6 wherein the means for calculating an 

2 opacity gradient for each of the voxels in the volume using the calculated revised values 

3 of opacity for all voxels further comprises: 

4 means for combining the revised values of opacity for the voxels in each cell to 

5 derive three orthogonal opacity gradient components for the voxel in the center of each 

6 cell; and 
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7 means for combining the three orthogonal opacity gradient components for the 

8 voxel in the center of each cell to derive an opacity gradient that is normal to an 

9 isosurface passing through the voxel in the center of each cell. 

1 18. The volume rendering apparatus of claim 14 wherein the means for calculating an 

2 opacity gradient for each of the voxels in the volume using the calculated revised values 

3 of opacity for all voxels further comprises: 

4 means for combining the revised values of opacity for the voxels in each cell to 

5 derive three orthogonal opacity gradient components for the voxel in the center of each 

6 cell; and 

7 means for combining the three orthogonal opacity gradient components for the 

8 voxel in the center of each cell to derive an opacity gradient that is normal to an 

9 isosurface passing through the voxel in the center of each cell. 

1 1 9. The volume rendering apparatus of claim 1 7 further comprising: 

2 means for calculating shading for the volume using the opacity gradient; and 

3 means for displaying the rendered volume on a display device. 

1 20. The volume rendering apparatus of claim 14 further comprising: 

2 means for calculating shading for the volume.using the opacity gradient; and 

3 means for displaying the rendered volume on a display device. 

1 21. Apparatus for rendering a volume of voxel data with shading and opacity 

2 dependent upon the direction of a source of light specified to be illuminating the volume 

3 for rendering, wherein each voxel comprises a value representative of a parameter at a 

4 location within the volume and each voxel has an initial value of opacity, the apparatus 

5 comprising: 

6 (1) means for calculating a revised value of opacity for a first voxel in the volume * 

7 by creating a cell of voxels surrounding the first voxel, wherein all voxels in the cell are 
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8 arranged into groups so that each voxel in each group of voxels within the cell are 

9 positioned one behind the other in a line parallel to the primary direction of light, the 
10 revised opacity calculation means further comprising: 

1 1 -P^^ «~*+;~,-r ~~,r:„~J , ~ -C - „ -T 1 ~1 A .L _ 

a jlxa^oaao iui o^tum^ cxiv/ ibvoou vamc ux uuauiiy ux lijlc vuaci uiuscsi lu 

12 the source of light in each group of voxels in the cell equal to its initial value of 

13 opacity; and 

14 (b) means for setting the revised value of opacity of other voxels in the 

1 5 groups of voxels in the cell equal to the revised value of opacity of an adjacent 

16 voxel in the same group of voxels that is closer to the source of light if the revised 

17 value of opacity of the closer voxel is equal to or greater than the initial value of 

1 8 opacity of the other voxel, and setting the revised value of opacity of the other 

19 voxel equal to the initial value of opacity of the other voxel if the revised value of 

20 opacity of the closer voxel is less than the initial value of opacity of the other 

21 voxel; 

22 (2) means for repeating the revised opacity calculations for all other voxels in the 

23 volume other than the first voxel; and 

24 (3) means for calculating a gradient of opacity for each of the voxels in the 

25 volume using the calculated revised values of opacity for all voxels. 

1 22. The volume rendering apparatus of claim 21 wherein the means for calculating an 

2 opacity gradient for each of the voxels in the volume using the calculated revised values 

3 of Opacity for all voxels further comprises: 

4 means for combining the revised values of opacity for the voxels in each cell to 

5 derive three orthogonal opacity gradient components for the voxel in the center of the last 

6 mentioned cell; and 

7 means for combining the three orthogonal opacity gradient components for the 

8 voxel in the center of each cell to derive an opacity gradient that is normal to an 

9 isosurface passing through the voxel in the center of the last mentioned cell. 
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3 


2 


23. 


The. volume rendering apparatus of claim 21 further comprising: 

means for calculating shading for the volume using the opacity gradient; and 

means for displaying the rendered volume on the display device. 


24. 


The volume rendering apparatus of claim 22 further comprising: 

means for calculating shading for the volume using the opacity gradient; and 

means for displaying the rendered volume on the display device. 


2 


1 25. Apparatus for rendering a volume of voxel data with shading and opacity 

2 dependent upon the direction of a source of light specified to be illuminating the volume 

3 for rendering, wherein each voxel comprises a value representative of a parameter at a 

4 location within the volume and each voxel has an initial value of opacity, the apparatus 
.5 comprising: 

6 (1) means for calculating a revised value of opacity for a first voxel in the volume 

7 by creating a cell of voxels surrounding the first voxel, wherein all voxels in the cell are 

8 arranged into groups so that each voxel in each group of voxels within the cell are 

9 positioned one behind the other in a line parallel to the primary direction of light, the 

10 revised opacity calculation means further comprising: 

1 1 (a) means for setting the revised value of opacity of the voxel closest to 

12 the source of light in each group of voxels in the cell equal to its initial value of 

1 3 opacity; 

14 (b) means for setting the revised value of opacity of other voxels in the 

1 5 groups of voxels in the cell equal to the revised value of opacity of an adjacent 

16 voxel in the same group of voxels that is closer to the source of light if the revised 

17 value of opacity of the closer voxel is equal to or greater than the initial value of 

1 8 opacity of the other voxel, and setting the revised value of opacity of the other 

1 9 voxel equal to the initial value of opacity of the other voxel if the revised value of 

20 opacity of the closer voxel is less than the initial value of opacity of the other 

21 voxel; 
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22 (2) means for repeating the revised opacity calculations for all other voxels in the 

23 volume other than the first voxel; 

24 (3) means for calculating three orthogonal opacity gradient components for the 

25 first voxel, and one of the orthogonal opacity gradient components is parallel to the 

26 primary direction of light, the opacity gradient calculation means further comprising : 

27 (a) means for combining the revised values of opacity for the voxels in 

28 the cell to derive the three orthogonal opacity gradient components for the first 

29 voxel; 

30 (b) means for combining the three orthogonal opacity gradient 

3 1 components for the first voxel to derive an opacity gradient that is normal to an 

32 isosurface passing through the first voxel; and 

33 (c) means for repeating the opacity gradient calculations for all other 

34 voxels in the volume other than the first voxel; and 

35 (4) means for calculating shading for the volume using the opacity gradient. 

1 26. The volume rendering apparatus of claim 12 further comprising: 

2 means for calculating shading for the volume using the opacity gradient; and 

3 means for displaying the rendered volume on the display device. 

1 27. A computer readable medium containing executable instructions for rendering on 

2 a display device a volume of voxel data with shading and opacity dependent upon the 

3 direction of a source of light specified to be illuminating the volume for rendering, 

4 wherein each voxel comprises a value representative of a parameter at a location within 

5 the volume and each voxel has an initial value of opacity, the executable program 

6 instructions comprising program instructions for: 

7 (1) calculating a revised value of opacity for a first voxel in the volume by 

8 creating a cell of voxels surrounding the first voxel, wherein all voxels in the cell are 

9 arranged into groups so that each voxel in each group of voxels within the cell are 
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10 positioned one behind the other in a line parallel to the primary direction of light, the 

1 1 revised opacity calculation means further comprising: 

12 (a) setting the revised value of opacity of the voxel closest to the source 

13 of light in each group of voxels in the cell equal to its initial value of opacity; and 

14 (b) setting the revised value of opacity of other voxels in the groups of 

15 voxels in the cell equal to the revised value of opacity of an adjacent voxel in the 

1 6 same group of voxels that is closer to the source of light if the revised value of 

17 opacity of the closer voxel is equal to or greater than the initial value of opacity of 

1 8 the other voxel, and setting the revised value of opacity of the other voxel equal to 

19 the initial value of opacity of the other voxel if the revised value of opacity of the 

20 closer voxel is less than the initial value of opacity of the other voxel; 

21 (2) repeating the revised opacity calculations for all other voxels in the volume 

22 other than the first voxel; 

23 (3) calculating a gradient of opacity for each of the voxels in the volume using 

24 the calculated revised values of opacity for all voxels; and 

25 (4) calculating shading for the volume using the opacity gradient. 

1 28. A computer readable medium containing executable instructions for rendering on 

2 a display device a volume of voxel data with shading and opacity, wherein each voxel 

3 comprises a value representative of a parameter at a location within the volume and each 

4 voxel has an initial value of opacity, the method comprising the steps of: 

5 calculating a revised value of opacity for each of the voxels in the volume, the 

6 revised value of opacity of a voxel being dependent upon its initial value of opacity and 

7 the revised value of opacity of voxels proximate thereto; 

8 calculating an opacity gradient for each of the voxels in the volume using the 

9 calculated revised values of opacity for all voxels; and 

10 calculating shading for the volume using the opacity gradient. 

$ $ * sfc $ 
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